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ABSTRACT 

The launch of Hinode satellite led to the discovery of rising plumes, dark 
in chromospheric lines, in quiescent prominences that propagate from large 
(~ 10 Mm) bubbles that form at the base of the prominences. These plumes 
present a very interesting opportunity to study Magnetohydrodynamic (MHD) 
phenomena in quiescent prominences, but obstacles still remain. One of the 
biggest issues is that of the magnetic field strength, which is not easily measur- 
able in prominences. In this paper we present a method that may be used to 
determine a prominence's plasma (3 when rising plumes are observed. Using the 
classic fluid dynamic solution for flow around a circular cylinder with an MHD 
correction, the compression of the prominence material can be estimated. This 
has been successfully confirmed through simulations; application to a prominence 
gave an estimate of the plasma (3 as /3 = 0.47±0.079 to 1.13±0.080 for the range 
7 = 1.4-1.7. Using this method it may be possible to estimate the plasma of 
observed prominences, therefore helping our understanding of a prominence's 
dynamics in terms of MHD phenomena. 

Subject headings: Magnetohydrodynamics (MHD), SumProminences 



Introduction 



Quiescent prominences are large clouds of relatively cool plasma supported against grav- 
ity by their magnetic field. It is known that th e temperature of quiescen t prominences is 



approximately 8000 K (iTandberg-Hanssenl Il995l ) and density ~ 10 cm ( iHirayamal Il986f ) 
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which is a decrease and increase of approximately two orders of magnitude respectively from 
the surrounding corona. Using this value for the temperature the pressure scale height can be 
calculated as A ~ 300 km, which is approximately 2 or ders of magnitude less th an the char- 
acteristic height of a quiescent prominence (~ 25 Mm iTandberg-Hanssenl Il995l ) . There are 
many reviews that describe the current understanding of the nature of quiescent prom inences 
(see, for example, lTandberg-Hanssenlll995l ; lLabrosse et al.ll2010t iMackay et al.ll2010l ). 



Obs ervations by the Sol ar Optical Telescope (SOT; iTsuneta et al.ll2008l ) on the Hinode 
satellite ( iKosugi at al.l 120071 ) have shown that on very small scal es qu iescent prom i nences 
are highly dynamic and unstable phenomena. iBerger at al.l ( 120081 ) and iBerger at al.l (120101 ) 
reported dark upflows that propagate from large (~ 10 Mm in size) bubbles that formed at 
the base of some quiescent prominences, through a height of approximately 10 Mm before 
dispersing into the backgroun d prominence mate rial. An example of these plumes is shown 
in Figure [TJ Observations by IBerger et al.l ( 120111 ) show that the bubble and plumes have a 
minimum temperature of 250, 000 K. The dark upflows were found to rise at constant velocity 
of ~ 20kms~ 1 . Often these plumes would separate from the cavity forming bubbles inside 
the prominence material. 
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Fig. 1. — Panel (a) shows a quiescent prominence observed on the 3-Oct-2007 at 04:16UT. 
Panel (b) shows a zoom-in of the region in the box in panel (a). The plume and the bright 
prominence plasma above it are clearly visible. 



Berger at al.l (120 lOf ) hypothesized that the observed upflows were caused by the mag- 
netic Rayleigh- Taylor instability, as a mixed mode perturbation with an interchange mode 
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(k perpendicular to B) and an undular mode (k parallel to B), in the high Atwood number 
limit [A = (p + — p_)/(p + + p_) — > 1], where fc is the wave number, p + is the density of 
the region ab ove the contact disco ntinuity and p_ is the density of the below the contact 
discontinuity. iRyutova et all ( 120101 ) described how the theoretically predicted growth rates 
and behaviour for the magnet ic Rayleigh- Taylor insta bility well match th e observations of 
quiescent prominence plumes. iHillier et al.l ( 120111 ) and iHillier et al.l (120121 ) showed how up- 
fiows could be created by the 3D mode of the magnetic Ray leigh- Ta ylor instability acting on 
the b oundary between the Kippenhahn-Schluter prominence model (jKippenhahn fc Schliiter 
19571 ) and a hot tube inserted into the it. The upflows created were found to be qualitatively 
similar to the observed plumes. 

Observations of th e magnetic fi eld of quiescent prominences show the field strength to be 
i n the range ~ 3- 30 G ( jLeroylll989l ). where the average for polar-crown prominences is ~ 5 G 
( lAnzer &: Heinzell 120071 ) . However, as it is very difficult to observe the prominence magnetic 
field directly, often other methods to determine the magnetic field of prominences have been 
employed. One common mechanism has been to extrapolate the photosp heric magnetic field 



to obtain values for the cor onal field where the prominence is observed (IDudik et al.l 12008 



Aulanier fc Demoulinl 20031) . Another common method is to analyse the oscilla tory pattern 



of filame nts (llsobe fc Tripathil l2006at iPinter et al.l l2008t iTripathi et al.l 120091 ) or filament 



threads ( ILin et al.l 120071 120091 ) to calculate the magne tic field through the application of 
wave theory. This is known as prominence seismology (lArregui et al.l 120121 ) . In this paper 
we present a new method to investigate the magnetic field strength relative to the gas pressure 
by determining the plasma /3 (ratio of gas pressure to magnetic pressure) of a prominence 
through the study of compression of the prominence material by the Rayleigh- Taylor plumes. 



2. The Observations Needed to Determine a Prominence's Plasma f3 

Before the methodology through which the prominence plasma /3 can be determined 
is presented, first we will quickly review the observational characteristics of the promi- 
nence plumes that are important for this analysis. The points we need to remember are 



f lBerger at al.ll2010h : 



The plumes have a phase in their evolution where they rise at an approximately con- 
stant velocity 

The plumes have an approximately circular (elliptical) head 

The intensity of the prominence material at the top of the plume is greater than the 
average prominence intensity. 
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Flows of material can be observed along the head of the plume 



To make it possible to use these observables, a number of assumptions are necessary. 
The first is in relation to the density of the prominence. The prominence intensity is known 
to be a result of scattering. Therefore the intensity is directly related to the column density 
of the prominence. This allows us to assume that the intensity can be used as a proxy for 
the density. Hence the bright top of the plume would imply that there has been an increase 
in the density through compression of the prominence. 

The next assumption is in relation to the modelling of the geometry of the problem. 
As the plume moves at an approximately constant velocity with an elliptical head, this 
can be used to simplify the system under study. The first step is to remove the stem of 
the plume, and just leave an ellipse in the prominence material. Then as this ellipse is 
moving at constant velocity through (what we assume to be) a constant medium, a shift 
in reference frame can make the situation of an ellipse, stationary in a constant flow. The 
effect of a non-stationary prominence will be investigated to some extent later. Now we must 
consider the 3D nature of the geometry under investigation. The observed plumes are known 



to be caused by the magnetic Rayleigh- Taylor instability ( iBerger at al.ll2010l ; iHillier et al. 



20121 ) It is also known that the 3D magnet ic Rayleigh- Taylor instability creates filamentary 
structures aligned with the magnetic field (jlsobe et al.ll2005l l2006bl ; IStone fc Gardinerll2007l ; 
Hillier et al.l 120121) or the bi-sector of the upper an lower magnetic fields in the case where 
shear is present (IStone fc Gardiner! 120071 ) . Therefore we can view the problem as being 
similar to a elliptical cylinder in a constant flow. Then a coordinate transform to cylindrical 
coordinates would effectively give a circle in a constant flow. This has now reduced to a 
problem that has been significantly studied in fluid dynamics. 

The transform described above is a conformal transform of the coordinate system. In 
an incompressible regime, such a conformal transform would automatically give the correct 
result, but in a compressible regime, where nonlinearity is important, the velocity will not 
be completely accurate. The smaller the eccentricity and the fast-mode Mach number, the 
more accurate the coordinate transform will be. 



3. Derivation of equations to determine a prominence's plasma (3 

Here we use a classic result from fluid dynamics to define a set of equations that can 
be used to interpret prominence observations and determine physical parameters relating to 
that prominence. First, we describe the thought process behind the method, followed by a 
derivation of the equations used to determine the plasma /3 of the prominence. 



-5 - 



The derivation presented here is analogous to that presented in Ivan Dyke! (119751 ) for a 
compressible flow around a circular cylinder in a steady state. Here the derivation is altered 
to include a magnetic field that is perpendicular to the plane in which the plasma is flowing. 
In this derivation we assume that the magnetic field is in the x - direction and the flow has ve- 
locity components in the y- and z- direction. All derivatives in the x- direction are assumed 
to be 0. The validity of these assumptions will be investigated through comparisons with 
simulations in section HI This analysis uses the equations of mass conservation, momentum 
and energy. 

The equation of mass conservation, for a steady planar two-dimensional flow is: 

where p, v and w are the density, the y - component of the velocity and the z - component 
of the velocity. The equation presented above is for a steady state solution, so there is no 
temporal variation of the density and the velocity field. All these values have been normalised 
by the free stream reference values for density and velocity Uoo, where the values with 
the subscript oo denote the values upstream of the plume head. This can then be rewritten 
to give: 

dv ^ dw v dp w dp 
dy dz pdy p dz 

If the flow is irrotational, then the momentum equations can be simply replaced by the 
statement that the vorticity is zero, that is: 

dw dv 
dy dz 

which is satisfied by the velocity potential where: 

d(p d<p 

The next step is to relate p to the velocity, so that Equation [2] can be expressed purely 
in velocity derivatives, that is to say purely in derivatives of 0. For a perfect gas under 
adiabatic contraction/expansion that has a magnetic field perpendicular to the plane of the 
flow, the MHD-Bernoulli's equation can be written as: 

\ (v 2 + w 2 ) + -^—^ + ^ = C(S) (5) 
2 v ' 7 - 1 p 4irp K ' w 

where C(S) is used to show that the value is constant along streamlines. Therefore, the 
values of the physical variables at each point in the flow can related to the free stream values 
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(the values of the flow at an infinite distance from the circular cylinder). 

^ ■ I(„W )f £ + ^k = ^_^ + V (6) 



7-1 ppoo 2 p inpoo 7 - 1 poo 2 47rp c 

where the subscript oo implies the free stream values that are used for normalisation. 

If we define the free stream speed of sound as C s = ^/7Poo/Poo, the free stream Alfven 
velocity as V a = B xoo / y/inpoo and denote the free stream Mach number as = U^/Cg 
then equation [6] becomes: 



1 p 1 . 2 2N 1 1 2 V? 

£ _i_ _m 2 fw 2 + w 2 ) H -— = h -M 2 + — 

7 - 1 p 2 p C 2 7 - 1 2 C 2 



s 



We can then define a relation between the Alfven velocity and the sound speed as follows: 



V\_ 2 



C 2 7/3oo 

where fi^ = Bl^/p^ is the free stream plasma (3. This then gives: 

1 P 1 w2 / 2 2x 5 2 2 1 1 ,„ 2 

P + -M^ 2 + w 2 ) + - = + -Ml + — - . (9) 



7-1P 2 ' P iPoo 7 - 1 2 00 7 /3 c 

As we have assumed that the flow is irrotational, we can assume that the flow is constant 
entropy. For a perfect gas (in normalised form) this can be expressed as: p/p 7 = 1. We can 
also use the identity that B 2 = p/j3', where P'Poo defines the plasma (3 of system. Using 
these identities, equation [9] becomes 

^ + + *t> - 1) + = 1 + (10) 

To tidy up the equations, we define: 

D( Vi z) = 1 + 1-^1 (11) 



/?' 7# 

(t: 

7/9« 



DC 



C=l + *^> (12) 



giving: 



^^nb-^^-^nh (13) 



For simplicity, we will now make the assumption that the plasma f3 is approximately 
uniform. In reality, we would have to solve for this using the induction equation, but for 
simplicity this is neglected in this analysis. This leads to the following equation: 

P 7 " 1 + ^^Ktf + w 2 -l) = l. (14) 
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Now it is possible to combine equation [TH with equation [2] to determine an equation for the 
velocity field. 



If we follow the derivation by Ivan Dykd ( 119751 ). by transforming to cylindrical coor- 
dinates (where the length-scale is normalised by the radius of the cylinder) we can then 
calculate the velocity potential (with all d/dx = 0): 



cos6> - Mj 
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where: 



(15) 



(16) 



which can be viewed as the fast-mode Mach number. Using v = V</>, we get 
1 



v r (r,6) 



cos 9 - MS 



13 1 3 1 5 1 

'l2r* + 2H ~ u7 6 
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+ 0(M* 4 ) 
1 



(17) 



+ 0{M, 



1 + — ) sin^ + M, 



13 1 11 11 
r H ^ I sin ( 

12 r 2 2r 4 12 r 6 
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3.1. Derivation of density distribution 

It is now necessary to derive an equation for the distribution of the density. We need 
to convert equation [14] to cylindrical coordinates: 

P 1 - 1 + 1 -^M*{vl + vl-l) = l (19) 

Rewriting Equation HH we find: 

P = f 1 - + vl - 1) J (20) 

This equation can now be used to determine the expected density distribution around a 
plume head. 

Figure [2] shows the density distribution around a circular cylinder. The high density 
region, created as a result of the compression of the plasma, is clearly visible at the top of 
the circular cylinder. This is for the case where the higher order terms were used, i.e. the 
distribution shown is from Equation [2TJ1 
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Fig. 2. — The density found for a compressible flow around a circular cylinder. Brighter 
colours imply higher densities. 



4. Comparison with Simulation Results 

As it is not possible to directly confirm the effectiveness of this method for determining 
the plasma /3 of a prominence, it is important to apply it to simu latio ns results to provid e 
some validation. We apply it to the results of iHillier et al.l ( 120111 ) and iHillier et al.l ( 120121 ). 
where a full description of the simulations are presented. In this section we apply Equation 
l2"Ulto solve the forward problem, i.e. we use the parameters of the simulation - including the 
plasma /3 - to see how well the density distribution around the plume head can be reproduced. 

The simulation used in this paper is of the 3D magnetic Rayleigh- Taylor instability 
in the Kippenhahn-Schiilter prominence model. To destabilise the Kippenhahn-Schiilter 
prominence model, a low density tube (density of 30% of the prominence) was placed inside 
the model. The boundary between this tube and the dense material was subjected to a 
random velocity perturbation allowing the magnetic Rayleigh- Taylor instability to develop. 
A 3D rendering of the initial conditions of this simulation are shown in Figure 

In this study, 3D conservative ideal MHD equations were used. Constant gravitational 
acceleration was assumed, but viscosity, diffusion, heat conduction and radiative cooling 
terms were neglected and an ideal gas was assumed. The equations were non-dimensionalised 
using the sound speed (C s = 13.2kms _1 ), the pressure scale height (A = C s /(yg) = 
R g T/(fig) = 6.1 x 10 7 cm), the density at the centre of the prominence (p(x = 0) = 10~ 12 



Fig. 3. — 3D rendering of the initial conditions of the simulation. The blue isosurface shows 
the density and the lines show the magnetic field. 



g cm 3 ) and the temperature (T = 10 4 K), giving a characteristic timescale of r = A/C s 



47 s. The ratio of specific h eats was set a s 7 = 1.05 and (3 
simulations can be found in Irlillier et al.l ( 2012 ). 



0.5. A full description of the 



4.1. Description of simulated plumes 

Panel (a) of Figure H] shows the density distribution at the centre of the the Kippenhahn- 
Schliiter model as the rising plumes, created by the magnetic Rayleigh- Taylor instability 
move upward through the dense prominence material. The plume in the box, more like a 
bubble in this image but really a low-density tube, is the plume we use for this comparison. 
Note that we didn't have to choose a plume that has detached from its tail, as we are only 
interested in the density at the head of the plume. 

Panel (b) of Figure H] shows the time-distance plot taken along the slit in Panel (a). The 
constant rise velocity of th e plumes, resulting from the force balance created at the head of 



the plume, is clearly shown iHillier et al.l ( 120121 ). The rise velocity of the topmost plume was 



found to have a Mach number of M = 0.28. We use the free-stream Mach number because 
if we transform the plume head to 0-velocity rest frame then the rise velocity transforms to 
a free-stream velocity of the dense material. 
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Fig. 4. — Panel (a) shows the rising plumes that are formed by the magnetic Rayleigh- Taylor 
instability in the Kippenhahn-Schluter model. Panel (b) shows the time-distance diagram 
taken along the dashed line in panel (a). This allows the Mach number of the rise velocity 
of the plume to be calculated as M = 0.28. 

It can be instantly noticed in Figure S] that the plume is not perfectly circular. The 
first step of the transform is to change the coordinate size so that the size of the elliptical 
bubble is that of a unit ellipse (i.e. the area of the ellipse is 7r). This is achieved by dividing 
the y and z axis by yfab, where a is the minor axis of the ellipse and b is the major axis 
of the ellipse. A second transform needs to be performed to transform the unit ellipse to a 
unit circle. This is achieved by multiplying the z axis by y/b/a and the y axis by \fafb. 
This process is shown in Figure The curved line marking the head of the plume is used in 
Figure [7]to compare the velocity distribution at the plume head to the expected distribution. 

4.2. Comparison of Simulated and Predicted Density Distribution and 



To apply Equation [2U1 various parameters have to be determined. The Mach number 
of the rising flow is shown in Figure H] to be 0.28. The density distribution needs to take in 
mass conservation, so the column density (integrating in the x-direction) needs to be used. 
We calculate the average density over x = 0- 2.5A, where x is the direction of the horizontal 
magnetic field in the simulation domain. In a similar method to determining the density 
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Fig. 5. — Calculation of the coordinate shift for the plume head, panel (a) shows the original 
size with the plume dimensions marked (a = 1.05A and b = 0.9A) and panel (b) shows the 
result of that transform. The arch at the top of the plume marks the size of the plume head. 

distribution, the average plasma (3 is calculated in the same way, giving (3 ~ 0.55. The ratio 
of specific heats (7) was taken as 7 = 1.05 for the simulations presented here. 

The next important step is to determine the background values. We determine the 
background density profile as the density profile at the same height at the dashed line along 
the left-hand edge of Figure |6j This is then removed from the distribution above the head of 
the plume to show the density increase due to the compression by the rising plume. Finally 
the background distribution, required to normalise the density distribution, is calculated 
from the average distribution (over x G [0,2.5A]) as p av ~ 0.67. One final point to note is 
that the background distribution has a weak upflow (Mach number M = 0.07) that needs to 
be included in the calculation to move the rising plume into a stationary reference frame with 
a free stream Mach number of = 0.21. It is important to note that because the frame 
of reference is changed, this free stream Mach number is considered as the Mach number of 
the flow toward a stationary plume head. 

Figure |6] shows the density profile of the simulation at the head of the plume, with 
contours from the predicted distribution to highlight the expected distribution. The high 
density region at the head of the plume and the rarefied regions at the sides of match well 
between the predicted and simulated densities. 



-2-10 1 2 

Y/b 



Fig. 6. — This shows the average density (averaging in the x-direction). The density increase 
at the head of the plume can be clearly seen. The white contours show the predicted 
distribution. The dashed black line shows the position of the slit used in Figure El The 
background distribution has been removed. 

Next we compare the velocity distribution at the plume head with the distribution 
predicted by Equations [IT] and [£E1 Figure [7J shows the predicted distributions (solid line) 
and the simulated distributions (dashed line) for the vertical and horizontal velocity. The 
horizontal scale uses the same scale as that of Figure but centers the zero position at the 
plume head. The velocity distributions are in general good agreement. It should be noted 
that the velocities are multiplied by a factor of p(x = 0)/ p av to give the correct magnitude 
of the velocity (i.e. the maximum value of W/Cs ~ 0.21 ). Note the good agreement in the 
velocity field even though the circular solution has been modified to apply to an elliptical 
body. 

Figure M shows the density distribution (minus the background distribution) at the top 
of a plume (solid line). The dashed line shows the expected density distribution when using 
the parameters discussed above to give a value of « 0.19. The plots demonstrate that 
the theoretically predicted density distribution and velocity field well describes those from 
the simulation. 

To investigate this further, we can fit the density distribution using Equation [2D] and 
solve for M* 2 . This, in turn, will make it possible to determine the value of the plasma (3. The 
method is explained in greater detail in the next section. The fit for the density distribution 



13 




Fig. 7. — Velocity distribution along the head of the plume (as marked in Figure [5] panel 
(a)). The solid line shows the simulated distribution and the dashed line shows the predicted 
distribution. 

gives M* 2 = 0.038, which gives a plasma /3 of f3 = 0.59 when solving for the unknown (i.e. 
/3) in Equation [16j This value is slightly higher than the input value, but sufficiently close 
to show that the method can estimate the plasma /3 using the compression of the plume. 
Figure [9] shows the chi-squared for the density distribution for models using different plasma 
values and the simulated distribution. The minimum of the x 2 distribution at a value 
that is close the value of plasma (3 in the simulation. 



5. Application to Observed Prominence Data 



In this section we apply the method presented above to determine the plasma (5 of the 
plume displayed in Figure [TJ In contrast to Section HI this is treated as an inverse problem, 
where using the intensity and rise velocity allows the plasma /3 of the prominence to be 
deduced. Due to the timescales of these compressions (~ 100 s) being shorter than the 
timescales for the cooling (~ 1000 s), we can assume to first-order approximation that the 
increased intensity comes through purely adiabatic effects. These adiabatic effects would be 
an increased column density and a temperature rise from the compression. Both of these 
can be calculated from the model described above. 



-14- 




Fig. 8. — Density profile along dotted line shown in Figure |6] minus the background density 
distribution (shown with the solid line) and the expected distribution using the parameters 
of the simulation explained in Section l4~2l (shown with dashed line). 



Figure [T] panel (a) shows a quiescent prominence seen on the NW solar limb on 2007 
October 03 observed by the SOT with the Ca II H filter at a cadence of 30 s. The time 
series of this observation was between 01:16UT and 04:59UT. This prominence presents 
many interesting dynamic features, for example the start of this observation (01:16UT) a 



l arge bubble has formed inside the prominence similar to those described in iBerger at al. 



( 120101 ). There are also a number of bright threads and downwardly propagating knots that 
occur durin g the duration of th e observations as well as the upwardly ejected plasma blobs 
presented in iHillier et all ( 120111 ). In this paper we are interested in the rising plume and, in 
particular, the bright rim that forms at the top of the plume. This is shown in panel (b) of 
Figure [TJ 



5.1. Calculation of the Plume Mach Number 

Before the plasma /? is calculated, it is necessary to calculate the Mach number of the 
rising plume. Following the method applied in Section I4.1[ we create a time-distance plot to 
show how the height of the plume changes with time. The dashed vertical line in Panel (a) 
of Figure [1] shows the position of the slit used. 
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Fig. 9. — x 2 distribution calculate for different values of plasma /3 giving different density 
distributions and the simulated distribution. 

Panel (a) of Figure [TOl shows the time-distance plot, where the bright, rising region is the 
rising lip of the plume. The cross marks show the top of the plume. This is calculated as the 
pixel with the greatest intensity change from the pixel below change multiplied by the pixel 
intensity. Using these pixels, the velocity of the plume and the velocity error were calculated 
to be Vpi ume = 12.3 ± 0.6 km s -1 . As this rise velocity is calculated from a time-distance 
diagram, it should be viewed as a lower limit on the velocity. The dashed line in the figure 
shows a velocity of ~ 12kms _1 . Throughout the rise of the plume, the rise velocity can be 
seen to be approximately constant. 

To calculate the Mach number of a flow it is necessary to know the sound speed 
of the ambient material. The sound speed of the dense prominence material is C s = 
[y(R/ fj^Tprom] 1 / 2 = llkms" 1 for 7 = 5/3, R = 8.3 x 10 7 ergK" 1 moU 1 , fi = 0.9 and 
Tprom — 8000 K. Using the velocity calculated above we find the Mach number of the flow to 
be Mqo = 1-12 ± 0.05 in the case where 7 = 5/3. 

Panel (b) of Figure fTUl shows the dimensions of the ellipse used. In this case a = 1.3 arcsec 
and b = 2.35 arcsec. These lengthscales can now be used to determine the normalising length 
scales for the system, so that the plume head becomes circular. Now the prominence plasma 
(3 can be determined. 
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Fig. 10. — Panel (a) shows time-distance diagram along the slit shown in Figured] panel (b). 
The cross marks show the position of the top of the plume. The rise velocity of the plume 
is found to be ~ 12kms _1 . Panel (b) shows the dimensions of the ellipse used to model the 
rising plume. 



5.2. Determining Plasma (3 of Observed Prominences 

Figure [11] shows the normalised intensity. Using the emissions for heights greater than 
4 in Figure E] a linear fit was used to de-trend the data (as was applied to the simulations). 
This intensity is used as a proxy for the column density assuming the emission is scattering 
dominate (oc p). In this paper we do not take into account the change in emission from the 
change in temperature. The fitting for the density is shown by the dotted line. This fit is 
performed using the IDL routine curvefit.pro. The fit solves for the value of M* 2 giving the 
values shown in Tabled] 

It is important to determine the error associated with this fit, so that the error in the 
plasma (3 can later be calculated. There are two important errors that are associated with 
this fit: the intensity error and the error between the observed curve and the model curve. 
These errors then determine the error for M* 2 . 

The intensity error is determined in a simple fashion. The fluctuations of intensity in 
the corona above the prominence are used. Here we assume that the fluctuations in the Ca 
II H line (which should not have a coronal signal) are due to stray light, and so we assume 
this is representative of the error in the intensity. For a region in the corona above the 
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prominence of size [200,100] pixel taken over 100 consecutive frames, the standard deviation 
of the intensity fluctuations is given as a cor = 1.33. As the fitting is normalised, the error 
should also be normalised. Performing this normalisation results in a normalised standard 
deviation of a cor = 0.028. This values is set as the intensity error of the prominence, i.e. 
0.028. It should be noted that this value is approximately the same as the standard 



a 



prom 



deviation of the fluctuations shown in the prominence in Figure [TT] for heights greater than 
4. 

The fitting error is determined as the standard deviation of the difference between the 
normalised intensity and the fit at each pixel. These errors are shown in Table [TJ The error 
bars equivalent to 2a prom and 2&f it are shown in Figure [TTJ The errors for M* 2 from the 
fitting routine are also shown in Table [TJ 




4 6 
Normalised length 



Fig. 11. — The graph shows a plot of intensity along the slit (solid line) and the fitting as- 
suming the intensity proportional to density (dashed line) for 7 = 1.65. Error bars equivalent 
to twice the standard deviation are presented for the observed intensity and the fit. 



Using the value of M* from the fitting of the intensity and the definition of M* as: 




7&o 

+ 2(7 - 1, 

It is possible to solve for the plasma /3 of the prominence, where (3 is given by: 

2(7 - 1) Ml 



7 



M 2 - M 2 



(21) 



(22) 
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The results for the plasma (3 for different values of 7 are given in Tabled] The range of values 
covers two orders of magnitude from j3 = 0.044 ± 0.078 for 7 = 1.05 to (3 = 2.666 ± 0.080 for 
7 = 2.0. However, the value for 7 in a prominence is likely to be in the range 7 = 1.4- 1.7, 
giving a range of plasma (3 of = 0.47 ± 0.079 to 1.13 ± 0.080. This is a range of only a 
factor of ~ 2.5. 



6. Discussion 

In this paper, using Equation [20] we found that for the range 7 = 1.4- 1.7 the observed 
prominence has a plasma [3 of (3 = 0.47 ± 0.079 to 1.13 ± 0.080. Using Equation [2171 and the 
analysis presented in this paper, it is possible to make some predictions on the nature of the 
density enhancement at the head of a plume: 

• The larger the value of M* (i.e. the larger (3 and are), the larger the increase 
in density at the head of the plume. This is a simple consequence of the greater 
perturbation from the incompressible state the system receives. 

• The larger the rising plume, the greater the thickness of the region at the head of the 
plume over which the density increases. This has important implications relating to 
which plumes this mechanism can be accurately applied. 

It is clear from Table [T] that 7 presents the biggest uncertainty relating to the value of 
the plasma (3 of the prominence. It may be expected that by looking at Equation [20] the fit 
may provide some constraints on the value of 7, but the of the fit is almost the same for 
each value of 7 therefore statistically it is not possible to do this. Separate methods that 
can constrain the value of 7 would work well in conjunction with this method to determine 
the plasma (3 of a prominence. 

The compressions that are under investigation here happen on timescales of approxi- 
mately 1 Mm/ 10km s -1 ~ 100s, which is shorter than the radiative cooling timescales of 
the prominence (~ 1000 s), so the compression should be well modelled as adiabatic. Cur- 
rently there is no model to describe accurately how the intensity of prominences changes 
under adiabatic compression. With this in mind, we will only use the column density, and 
ignore the effect of heating. We can at least be confident that the heating will not move the 
prominence plasma out of the observable temperature range, as the temperature increase is 
given as T' = p' 7 " 1 ~ p' 2 / 3 when 7 = 5/3. Therefore, an increase of the density of ~ 1.4 
(which is approximately what we find in this paper), would increase the temperature from 
8000 K to 10 4 K which is still in the emission range for Ca II H ~ 5000- 10 4 K so would not 
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Table 1: Results for Ml and plasma (with errors) for various values of 7 



7 


2 

X 


or fit 


Mt 









w) 


1.05 


1 


.410 


0.61 


0.027 


0.04 





,078 


1.1 


1 


.407 


0.62 


0.027 


0.09 





,079 


1.15 


1 


.404 


0.63 


0.027 


0.14 





,079 


1.2 


1 


.401 


0.63 


0.028 


0.20 





,079 


1.25 


1 


.400 


0.64 


0.028 


0.26 





,079 


1.3 


1 


.394 


0.65 


0.029 


0.32 





,079 


1.35 


1 


.390 


0.65 


0.029 


0.39 





,079 


1.4 


1 


.386 


0.66 


0.030 


0.47 





,079 


1.45 


1 


.382 


0.67 


0.030 


0.56 





,080 


1.5 


1 


.377 


0.67 


0.030 


0.64 





,080 


1.55 


1 


.372 


0.68 


0.036 


0.74 





,080 


1.6 


1 


.367 


0.68 


0.031 


0.86 





,090 


1.65 


1 


.362 


0.69 


0.031 


0.99 





,080 


1.7 


1 


.357 


0.70 


0.031 


1.13 





,080 


1.75 


1 


.351 


0.71 


0.032 


1.30 





,080 


1.8 


1 


.345 


0.71 


0.032 


1.49 





,080 


1.85 


1 


.339 


0.72 


0.032 


1.72 





,080 


1.9 


1 


.333 


0.73 


0.033 


1.98 





,080 


1.95 


1 


.327 


0.74 


0.033 


2.29 





,080 


2.0 


1 


.321 


0.75 


0.033 


2.67 





,080 
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move the emission out of the passband. A model of the intensity change due to the increase 
in temperature under these conditions would increase the accuracy of the calculation from 
plasma (3. 

The value for plasma (3 found here has a strong dependence on the hydrodynamic Mach 
number as the values found were not small enough to make it unimportant. Therefore, if 
downflows of prominence material interact with the head of the plume, they would increase 
the Moo, which would reduce the value of the plasma (3. In section loTTI we assume a static 
prominence, but if we use a prominence ba ckground w here the material is falling at sound- 



speed (based on the downflows observed by [Chad 12010), we can calculate a new value for the 



plasma (3 of (3 ~ 0.14 for 7 = 1.65. Figure [10] does not show any obvious downflows inter- 
acting with the head of the plume, so the assumption of ~ 11 kms" 1 downflows is probably 
a huge over estimate, but this calculation is shown to highlight the potential importance of 
downflows when using this method to calculate the plasma (3. Therefore, when downflows 
are not included in a calculation, as is the case in this paper, then the plasma (3 value should 
always be viewed as an upper limit. 

Another point that should be noted is the assumption of an irrotational flow as used in 
Equation [3j The most important point is that the plume is rising through the prominence 
material, pushing it out of the way as it rises. In a compressible medium with finite Reynolds 
number, the downstream flow (flow behind the plume head) would develop turbulence and 
a wake, breaking the irrotational assumption. However, the compression that is being mod- 
elled is at the head of the plume, where the flows are much less complex. There are also 



observations that suggest prominences are weakly turbulent (ILeonardis et al.l 120121 ). which 
would imply that to a some extent the irrotational assumption is also broken upstream of 
the plume. But it can be expected that to first order the properties of the mean flow are 
not affected by the turbulence, but the turbulence would be modified by the mean flow po- 
tentially changing the characteristics of the turbulence downstream of the plume head, i.e. 
turbulence distortion. To summarise, the assumption used is broken downstream, but holds 
to first order in the area we are interested in. Therefore, even though we couldn't use this 
to describe all the flows associated with the prominence plumes, it does apply for the area 
that is of interest allowing the plasma [3 to be determined. 

One interesting point to note is that we have equations that define the velocity around 
the plume head, as analysed in Figured Therefore, analysis of the velocities at the observed 



plume head (either thro ugh Doppler shift analysis or the cork-tracking used in lBerger at al. 



(120101 ) and IChad (120101 ) will allow the projection of the plume on the plane-of-sky to be 
determined. This would be important as it would remove some uncertainty in the width 
of the plume head, making the parameters of the fitted ellipse more accurate. With this 



-21 - 



information, it should be possible to determine the prominence's plasma f3 with greater 
confidence and also, through application of knowledge of the 3D magnetic Rayleigh- Taylor 
instability, estimate the direction of the magnetic field in the prominence. We plan to present 
this work in a future paper. 



It should be noted that ICargill et al.l ( 119961 ) studied a similar geometry through nu- 
merical simulations to investigate the movement of a magnetic cloud through a magnetised 
atmosphere in relation to CME propagation. Therefore, the method presented in this paper 
may be applicable to more MHD phenomena than just prominences, an investigation into 
what phenomenon this can be applied to may open up some very rewarding research areas. 
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